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1 Introduction 

A multitude of absorption lines seen shortward of QSO Ly-a emission, that cannot be traced to heavy 
element absorption systems, are assumed to be Ly-a lines arising in intervening clouds. Studies of these Ly-a 
clouds, typically done at 1 A or lower resolution, have shown N(HI) ~ 10 13 — 10 17 cm~ 2 and b ~ 35 km/s. 
Sargent et al 1980, on the basis of a flat pair velocity correlation function (PVCF), argued that these clouds 
are intergalactic. But Crotts 1989 showed that the strong Ly-a lines are spatially clustered. High resolution 
studies of Webb 1987 and Rauch et al 1992 also report some evidence for weak clustering, but overall such 
high resolution studies have been rare. Here we report a study of the Ly-a forest of Q1331+170 over 
z a bs — 1.60 — 2.19 based on i8 km/s resolution data at S/N ~ 15 , with metal-line deblending incorporated. 

2 Observations and Sample Generation 

The study of Q1331+170 ( z em = 2.08) made by York et al 1992 consisted of : (a) 18 km/s, mean S/N ~ 
15 Kitt Peak Echelle Spectra (over 3170-3970 A ) (b) 1 A, S/N ~ 40 MMT spectra (3900-9400 A) (c) limited 
35 km/s resolution MMT scans over 6400-6820 A. Heavy element lines, including those in a previously 
known damped Ly-a system at z 0 j s = 1.7765, were identified and analyzed (York et al 1992). On the basis 
of analysis of heavy element lines from known systems longward of the emission Ly-a, we derive contributions 
of heavy element lines to lines in the Ly-a forest using profile fitting techniques. The remaining contribution 
is then the deblended Ly-a lines. The sample of lines thus derived is further searched for previously unknown 
metal line systems. Two samples are derived from the analysis of the resultant list of pure Ly-a lines : one 
consisting of single component profile fits (SI) and the other consisting of the minimum required number of 
multiple components (S2). Unidentified lines slightly longward of the presumed QSO redshift ( z em = 2.08) 
are also included. Sample SI consists of 83 and sample S2 of 124 lines between z„b , = 1.60 — 2.19. Compared 
to a sample which would have ignored the metal-blended Ly-a lines completely, our 'procedure of metal 
line deblending adds ~ 20% lines to each of the samples SI and S2, and ~ 23% and ~ 29% lines to the 
sub-samples of strongest lines from SI and S2, respectively. 

3 Equivalent widths, Column densities, Doppler parameters 

The highest S/N is achieved between ~ 3500 — 3880 A resulting in rest equivalent width sensitivity of 
20 mA (4-5 a) in the best parts of the spectrum. At this sensitivity, most Ly-a lines are found to be weak. 
Distributions of HI column densities for samples SI, S2 are shown in Fig. 1. Only a very small fraction 
of lines. from sample S2 possess N(HI) > 10 14 cm -2 . This is different from the power law distribution 
(~ N(H I)~ 17 , log N (HI) > 13), which is common in the literature. (See for example, Carswell et al 1984.) 
According to this distribution, one would have expected about 9 lines with N(HI) > 10 14 5 in sample S2, 
whereas we see none. This is probably a result of the high resolution of our observations. High resolution 
studies of Ly-a lines in more sightlines would be needed to verify the generality of this. The b values range 
between 10 and 40 km/s for most of the lines, with a mean of ~ 27 km/s. No obvious correlation is apparent 
between N(HI) and b. Lines within ± 10,000 km/s of the presumed QSO emission redshift, z em = 2.08, are 
predominantly weaker than 100 mA, while lines ‘outside’ 10,000 km/s of the QSO are relatively stronger. 
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4 Velocity Correlations 

Fig. 2 shows the inner 8000 km/s of the histogram of ‘corrected’ number of line pair separations for all 
lines in sample S2 with W rett > 30mA, for go = 0.5. Here, ‘corrected’ number means the number of pairs 
after correcting for the finite redshift range of the data by using a ramp-shaped function, as in Sargent et al 
1980. Shown overlaid on the histogram are the mean and the ±2? levels expected in samples of randomly 
distributed Ly-a lines, computed as average of 100 simulated linelists, each having the observed number of 
lines with the observed set of equivalent widths, but distributed randomly between 2 0 » , = 1.6079 and 2.1914. 
An excess in number of velocity pairs out to 100 km/s is apparent in the Ly-or lines toward Q1331+170. This 
result seems to differ from that of Rauch et al 1992, who, with observational setup almost identical to ours 
and data at 23 km/s resolution, did not report such a clustering. This could be a result of different methods 
of analyzing the data. However they do report clustering among the narrow Ly-or lines. We also find a 
marginal excess in the number of line pairs over v <~ 500 km/s in the weakest (30 < Wrest < 100mA) 
lines, and the two findings could be related. An excess in number of velocity pairs out to ~ 200 km/s is also 
seen in the (admittedly small) sample of the strongest (W rett > 190mA) lines. This excess comes entirely 
from a strong clump of Ly-a lines at 2 a j, = 1.9637. 

5 Conclusions 

We have pruned the Ly-a forest of Q1331+170 with metal line deblending and find that our sample of 
Ly-a clouds have low N(HI) and show an excess in the number of tine pair separations on scales of <~ 200 
km/s. We have searched for CIV doublets and OI A1302-SUI A 1304 pairs, since these have Av « 500 km/s. 
We found only one, which has been taken out of the samples. The implied excess in the PVCF could be 
either due to clustering of the Ly-a clouds or due to complex structure in the clouds. The former conclusion 
is not expected in an intergalactic interpretation of the Ly-a clouds. This effect will have to be tested by 
more high resolution, high S/N studies of the Ly-a forest in other QSO’s. 

This work is part of a more detailed paper by Kulkarni et ml 1992. 
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Fig. 1 Distribution of N(HI): Samptos SI. S2 
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